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*^We describe the results of our recent work on the determination of the value of the parameter 
and of the Q 2 -dependence of the Gross-Llewellyn Smith (GLS) sum rule from the experi- 

■ Mb 

mental data of the CCFR collaboration on neutrino-nucleon deep-inelastic scattering, using the 
Jacobi polynomials QCD analysis. The obtained results are compared with the information, 
available in the literature, information on the previous experimental measurements of the GLS 
sum rule. 
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1. The detailed experimental and theoretical studies of the deep-inelastic scattering pro- 
cesses provide the important information on the applicability of perturbative principles for de- 
scribing the observed Q 2 -dependence of the nucleon structure functions (SF) Fi(x, Q 2 ) (i = 2, 3) 
and of the related moments 

M n (Q 2 ) = f 1 x^F^x.Q^dx (1) 
Jo 

within the framework of QCD. The non-singlet (NS) SF xF 3 (x, Q 2 ) = (xF 3 vN + xF£ N ) /2, which 
characterizes the difference of quark and antiquark distributions, can be measured in the deep- 
inelastic processes with charged electroweak currents. The most precise experimental data for 
this quantity were recently obtained by the CCFR group at the Fermilab Tevatron [|TJ. The 
comparison of these data with the perturbative QCD predictions for xF 3 was originally made 
with the help of the computer program developed in Ref. [[J, based on the solution of the 
Altarelli-Parisi equation. The fits were made for various Q 2 cuts of the data. In particular, 
fitting the data at Q 2 > 10 GeV 2 , the CCFR collaboration obtained the following value of the 



parameter 



(4) 

- Hill - 

(4) 



A^L = 171 ± 32(stat) ± 54(syst) MeV . (2) 

This value turns out to be almost non-sensitive to the variation of the Q 2 cuts, imposed for 
allowing one to neglect the effects of the high-twist (HT) contributions at low energies. 

Another important result, obtained by the CCFR collaboration, is the accurate measure- 
ment of the Gross-Llewellyn Smith sum rule 

GLS(Q 2 ) = - ^ ^tldx (3) 

2 jo x 

at the scale Q 2 = 3 GeV 2 0: 

GLS(Q 2 = 3 GeV 2 ) = 2.50 ± 0M8(stat.) ± 0.078(syst.). (4) 

It was already shown [[| that this value of the GLS sum rule is consistent with the QCD 
predictions, provided one takes into account not only the perturbative QCD corrections [||. Q 
to the quark-parton prediction GLSa s — 3, but the non-perturbative three-point function QCD 
sum rules estimates of the HT contributions as well. However, the interesting question of 



the possibility of extracting the Q 2 -dependence of the GLS sum rule from the CCFR data 
remained non-studied. 

In recent work we analysed this problem with the method of the SF reconstruction over 
their Mellin moments, which is based on the following expansion of the SF over the Jacobi 



polynomials P|- [pH| : 



Nmax = l2 U 

xF 3 N -*= 12 (x,Q 2 ) = x a (l-xf £ e^(x)E4 n) («,/5)M^ 2 (Q 2 ), (5) 

71=0 j=0 

where c^(a, 0) are the coefficients that are expressed through T-functions and the values of the 
Darameters a. 3. namelv a = 0.12 and 3 = 2.0. were determined in Ref. UTOl. The Q 2 evolution 



of the moments Mj^| can be determined from the solution of the corresponding renormalization- 
group equation expressed in the form presented in Ref.|l2|. The basic expansion parameter is of 
course the QCD coupling constant a s , which can be expressed through the QCD scale parameter 
in the standard way: a s (Q 2 )/4vr = ln(Q 2 /Af^) - fa lnln(Q 2 /Af^)//3 3 ln 2 (Q 2 /A§^) 
where (3 = 11 — 2/3/, /3 X = 102 — 38/3/. The relation of Eq. (5), supplemented by the 
corresponding solution of the renormalizat ion-group equation for M^ 2 , forms the basis of the 
computer program created by the authors of Ref. |TTJ. It was previously tested and used by 
the members of the BCDMS collaboration in the course of a detailed QCD analysis of the 
experimental data for F%{x, Q 2 ) SF of the deep-inelastic muon-nucleon scattering [13|. We were 
using in our studies also this program, thus building the bridge between the determination of 
AgL from F 2 (x, Q 2 ) and xF 3 (x, Q 2 ) SFs. 

2. In accordance with the original NS fit of the CCFR collaboration ]T|, ]!|, we have chosen 
the parametrization of the parton distributions at fixed momentum transfer in the simplest 
form 

xF 3 (x,Q 2 ) = A{Ql)x b ^\l - xf^\ (6) 

The constants A(Qq), b(Ql) and c(Qq) in Eq. (^) and the QCD scale parameter A^- were 
considered as free parameters, which were determined for concrete values of Q^. In order to 
avoid the influence of the HT effects and the TM corrections, following the original CCFR 
analysis we used the experimental points of the concrete CCFR data in the plane (x, Q 2 ) with 
0.015 < x < 0.65 and 10 GeV 2 < Q 2 < 501 GeV 2 . Note that we restricted ourselves by taking 
/ = 4 throughout the whole work. Moreover, we did not take into account any threshold effects 
in the process of our analysis. 

Using Eq. (5) we reconstructed the theoretical expression for xF 3 Nmax=12 (x, Q 2 , A, b, c, A) in 
all experimental points (x exp , Q 2 exp ) (for the detailed description see Ref. ||). The determination 
of the free parameters of the fit (namely A, b, c, A) was made by minimization of x 2 by the 
MINUIT program, which also automatically calculated their statistical errors. The numerical 
value of the GLS sum rule at different values of the reference scale Qq was determined by 
substituting the concrete values of the parameters A(Q 2 ) ), b(Ql) and c(Qq) into Eq. (7) and 
calculating its first moment, which determines the expression for the GLS sum rule. The 
statistical errors for the sum rule were calculated from the statystical errors of the parameters 
A(Q 2 ), b(Q 2 ) and c(Q 2 ). 

The results of the concrete calculations, made for various Qq points, are presented in 
Table 1. Figure 1 demonstrates our results for the GLS sum rule obtained in the process of 
the next-to-leading order (NLO) fit. It is worth emphasizing that putting n=l in the NLO 
expression for the moments M^ S {Q 2 ) one can reconstract the LO expression for the GLS sum 
rule only, namely GLSlo{Q 2 ) — 3[1 — a s /7r]. Threfore, in Fig. 1 we compare our result of the 
NLO fit with the LO perturbative expression of the GLS sum rule. 
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Table 1. The results of the LO and NLO QCD fit of the CCFR xF 3 
SF data for / = 4, Q 2 > 10 GeV 2 , N max = 12 with the corresponding 
statistical errors. The symbol Xd.f. * s f° r the x 2 parameter normalized to 
the number of degrees of freedom d.f. 




Fig.l. The comparison of the result of the NLO fit of the Q 2 evolution of the GLS sum rule 
with the statistical error bars taken into account with the LO perurbative QCD prediction. 



The value of the parameter A^™ in the LO perturbative expression for the GLS sum rule 

Ad O 

in Fig. 1 was taken in accordance with the results of our analysis of the CCFR data for the 
SF xF 3 at the reference point Ql = 3 GeV 2 (see Table 1). Notice, that the points presented in 
Fig. 1 are strongly correlated. The explanation is very simple: they were all obtained from the 
whole set of data. 

3. Taking into account our estimates of the statistical uncertainties and the estimate, 
determined by the CCFR group, of the systematic uncertainty H], we obained the following 
value of the GLS sum rule at the scale Ql = 3 GeV 2 ||: 



GLS(Ql = 3 GeV 2 ) = 2.446 ± 0.081(stat) ± Q.Q78(syst). 



(7) 



This in the agreement with the result (4) obtained by the CCFR group. The smaller statistical 
error of the CCFR result of Eq. (4) comes from their more refined analysis of this type of 
uncertainties. The results of our NLO determination of the GLS sum rule (see Table 1) do not 
contradict either the previous, less accurate, measurements of this sum rule at different values 
of Q 2 , as presented in Table 2. 
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10 
1-10 
10-20 

3 

3 

0.5-10 
3 
3 
3 



Result 



3.20 ± 0.50 
2.56 ± 0.41 ± 0.10 
2.89 ± 0.33 ± 0.23 
3.13 ± 0.48 ± 0.28 
2.83 ± 0.15 ± 0.10 
2.70 ± 0.40 
3.10 ± 0.60 
2.78 ± 0.08 ± 0.13 
2.50 ± 0.02 ± 0.08 
2.45 ± 0.08 ± 0.08 



Table 2. The summary of various determinations of the GLS 
sum rule with the corresponding statistical and systematical 
uncertainties. 



Note, that the results of the CHARM collaboration |L5| were obtained after integrating xF 3 SF 
in the interval 0.0075 < x < 1. We expect, that after interpolation of the data in the region 
of small x and integrating them in the whole interval < x < 1 the corresponding CHARM 
results for the GLS sum rule will approach the quark-parton prediction GLSas = 3 and will be 
even less accurate than the claimed result |T5j presented in Table 2. 

The results for A^- obtained from the NLO fit of the CCFR data, using the Jacobi poly- 
nomial expansion (see Table T), are in exact agreement with the outcome of the fit of the 



BCDMS data for the F 2 (x, Q 2 ) SF with the help of the same computer program [[0|, namely 
A^L = 230 ± 20(stat) ± 60(syst) MeV. 



The result, used at Fig. 1, namely A™. = 213 ± 31(stat) MeV , which was obtained 



MS 



using the expressions for the Mellin moments M^ s , is somewhat larger than the result of 
Eq. (2) obtained by the CCFR group with the method based on the solution of the Altarelli- 
Parisi equation. A similar feature was previously observed in the process of the analogous 
fits of the xFs less precise data obtained at Protvino [Ell]: the Altarelli-Parisi method gave 



A^L = 170 ± QO(stat) ± 120(syst) MeV (compare with Eq. (2)), while the fit over the Mellin 
moments resulted in the value A-^ = 230 ± 40(stat) ± 100(syst) MeV, which should be 
compared with the results of our fit (see Table 1). This slight discrepancy between the central 
values of the outcomes of different fits of the same data might be due to the intrinsic features 
of the two different ways the perturbative QCD corrections are taken into account. One can 
hope that this difference will be minimized after incorporating higher-order perturbative QCD 
effects into both methods. 

4. The Q 2 dependence of the GLS sum rule (see Fig. 1), extracted by us from the CCFR 
data, does not contradict the previous measurements of the Q 2 dependence of this sum rule, 



made by the BEBC-Gargamelle collaboration |T6| (see Fig. 2) and WA25 collaboration [18 
(see Fig. 3). 




Fig.2. Data on the GLS sum rule from the combined BEBC narrow-band neon and GGM-PS 
freon neutrino/antineutrino experiments [lfj. Errors shown are statistical only. 
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Fig.3. The Q 2 -dependence of the GLS sum rule, extracted by WA25 collaboration [|I 



However, the larger uncertainties of these data did not allow one to reveal the characteristic 
behaviour of the experimental results in the low-energy region, clearly seen in the analysis 
of the CCFR data || j8j. Indeed, the results of Eqs. (4) and (10) lie much lower than the 
theoretical predictions of the pure perturbative QCD. This feature demonstrates the importance 
of taking account of the HT contributions in the theoretical expression for the GLS sum rule. 
Their general structure is known from the results of Ref.|E§. The corresponding numerical 
calculations of these terms were made in Ref. [^| and more recently in Ref. ||23|| , using the 
same three-point function QCD sum rules technique. Combining all available information about 
the GLS sum rule we can write the following theoretical expression: 



GLSqcd{Q a 



I- a- (4.583 - 0.333/)a 2 - (41.441 - 8.02/ + 0.177/ 2 )a 3 - 



8 ((O)) 
27 Q 2 



where a = a s /7r is the coupling constant in the MS scheme. The NLO and NNLO perturbative 
QCD corrections were calculated in Refs. || and respectively. Note, that the perturbative 
expression for the GLS sum rule has one interesting feature, namely it is related to the per- 
turbative expression for the e + e~ —>■ hadrons D-function, calculated at the NLO and NNLO 
levels in Refs. [^4] and |25| (see also |26[), by the non-trivial generalization of the quark-parton 



connection of Ref. |27j • This generalization, discovered in Ref. |28| , has the perturbative cor- 
rections starting from the NLO level. They are proportional to the two-loop QCD /5-function. 
The theoretical consequences of the appearence of this factor are not yet clear. However, even 
at the current level of understanding it is possible to conclude, that the results of the work 
|28| provide the strongest argument in favour of the correctness of the results of the NLO and 
NNLO calculations of the GLS sum rule and the e + e~-annihilation R-ratio. 



Let us return to the discussions of the effects of the HT-contributions in Eq. (12). The 
original calculation Jj£| of the matrix element ((O)) gave the following estimate ((O)) = 0.33 ± 
0.16 GeV 2 . The QCD prediction of Eq. (12) with this value of the HT term was used in Ref. 
0| for the extraction of the value of from the experimental result of Ref. ||. We will 
not discuss here all details of the work of Ref. H, but present only final outcomes of the NLO 



analysis in the MS scheme with HT terms 



MS 

and without HT terms 



A§L = 318 ± 23(stat) ± 99(syst) ± 62(twist) MeV (9) 



AgL = 435 ± 20(stat) ± 87(syst) MeV . (10) 

It can be seen that the HT terms are decreasing the difference of the extracted values of the 
parameter from the results of other NLO fits, say from our results of Table 1. Even 
better agreement can be obtained after taking into account NNLO corrections in the GLS sum 
rule (see Eq. (12)), scheme-dependence ambiguities (see Ref. ||) or the new estimates of 
the HT-contribution, namely ((O)) = 0.53 ± 0.04 GeV 2 P3"|, which however have surprisingly 



small error bars of over 10% (it is known that the typical uncertainties of different three-point 
function QCD sum predictions lie within error bars of at least 30% ). 

Since at the scale Q 2 = 3 GeV 2 the result of our extraction of the GLS sum rule value || is 
in agreement with the original result of the CCFR group || (compare Eq. (4) with Eq. (10)), 
the conclusions of Ref. j|] remain valid in our case also. Moreover, we consider the deviation 
of the Q 2 dependence of the GLS sum rules results that we observed in the low-energy region 
from the prediction of perturbative QCD (see Fig. 1) as an indication of the necessity for a 
detailed study of the HT effects in the region of Q 2 < 10 GeV 2 . This conclusion joins the 
results of the quantitative analysis [Q, J29] of the effects of the HT contributions to the GLS 



sum rule |28 and the Bjorken polarized sum rules |30| correspondingly (see also Ref. [p3| ) and 



support the necessity of taking into considerations of these effects in the detailed description 
of the Q 2 dependence of the deep-inelastic scattering sum rules in the low energy region. 

In the high-energy region Q 2 > 10 GeV 2 the Q 2 behaviour of the GLS sum rule, obtained by 
us and depicted in Fig. 1, is in qualitative agreement with the perturbative QCD expectations. 
However, at the quantitative level there are indications of the existence of the deviation between 
theoretical predictions and the results of our analysis. This phenomenon might be related to 
the necessity for detailed studies of the effects of the NNLO corrections to the NS anomalous 
dimensions (which are known at present only for even moments |Jl[]) and of the NNLO coeffi- 
cients of the NS moments of the xF% SF ||32|| . Another important task is an improvement of the 



understanding of the behaviour of the xF 3 SF in the region of small x. We hope that a possible 
future analysis will allow one to study the Q 2 dependence of the GLS sum rule in more detail 
and to understand the status of the non-standard theoretical explanation of the behaviour of 
the GLS sum rule observed by us at moderate Q 2 p3 . 
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